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Abstract 

We investigate spatial behaviors of the CP-violating angle 9 by solving the equa- 
tion of motion of the two-Higgs-doublet model in the presence of a small explicit CP 
breaking 5. The moduli of the two Higgs scalars are fixed to be the kink shape with 
a common width. In addition to solutions 8 ~ 0(5) in all the region, we find several 
sets of two solutions of opposite signs, whose magnitudes become as large as O(l) 
around the surface of the bubble wall while the CP violation in the broken phase 
limit is of 0(8). Such set of solutions not only yield sufficient amount of the chiral 
charge flux, but also avoid the cancellation in the net baryon number because of a 
large discrepancy in their energy densities driven nonperturbatively by the small 5. 
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1 Introduction 



The electroweak baryogenesis |l[] strongly depends on the CP- violating angle 9(z) created 
at the first-order phase transition via bubble nucleation, where z is the spatial coordinate 
perpendicular to the 'planar' bubble wall. As is well known, one of difficulties of the 
electroweak baryogenesis in the minimal standard model may be that the explicit CP- 
breaking parameters of the Kobayashi-Maskawa scheme are too small to generate the 
observed baryon asymmetry of the universe. 

In a previous article |2|][], we examined the behaviors of 9{z) that is spontaneously 
generated in the two-Higgs-doublet model, by assuming that the moduli of the two neutral 
scalars, Pi(z)(i = 1,2), have the kink shape of a common width. One of the solutions we 
presented has a remarkable spatial dependence such that 9(z) becomes as large as 0(1) 
around the surface of the bubble wall while it completely vanishes deep in the broken and 
symmetric phase regions. Such a solution may be of much importance in the electroweak 
baryogenesis since, as we showed there, it does yield a large amount of the chiral charge 
flux through the wall surface at the phase transition. 

In this article we examine 9(z) in the two-Higgs-doublet model in the presence of an 
explicit CP breaking 5 at the transition temperature Tc, also by fixing the moduli to be 
the common kink shape. The magnitude of the breaking parameter S may not be largely 
different from those of the Kobayashi-Maskawa scheme even when finite-temperature cor- 
rections are taken into account. Then a naive guess would be that 9(z) remains of 0(5) 
in all the region between the broken and symmetric phase limits. Actually we give such 
solutions obtained analytically. On the other hand, we find several solutions whose 9(z) , s 
become as large as 0(1) around the wall surface, as important as the one in the case of 
spontaneous CP violation(<5 = 0). 

As pointed out by Comelli eta/.Q, the explicit CP breaking may be necessary to 
avoid the complete cancellation in the net baryon number expected from the symmetry 
of the solution 9(z) < — > —9(z) in the case of 5 = 0. In the presence of 5 ^ 0, the energy 
density of the bubble with 9 + (z) close to 9 s=0 (z) and that with 9~(z) close to —9 5=0 (z) 
no more degenerate. We give an estimate that the relative enhancement factor due to the 
energy difference between the two kinds of bubbles could be as large as O(10) even for 
5 ~ O(10~ 3 ). Such a large relative enhancement factor would favor the formation of one 
of the two kinds of bubbles and guarantee the baryon asymmetry of the universe. 

In Section 2 we introduce the breaking parameter 5 into the standard two-Higgs- 
doublet potential, and give the equation of motion for 9(z). In Section 3, we discuss the 

1 Ref.[|) is referred to as I. 
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boundary conditions to be satisfied by 9(z). In Section 4 we show examples of 9(z) ~ 
0(5) obtained analytically Several solutions of 9(z) ~ 0(1) around the bubble wall are 
presented in Section 5. In Section 6, the energy difference and the relative enhancement 
factor are estimated. Section 7 is devoted to concluding remarks. 

2 Explicit CP Breaking and Equation for 6 

In order to clarify essential roles played by 5, we examine the problem under the following 
simplified conditions: 

(1) One breaking parameter 5 is introduced into the standard two-Higgs-doublet poten- 
tial as (m|e~* <5 ${$2 + h.c.) while A5, A6, A 7 G R is assumed. 

(2) The magnitude of 5 at Tc is small enough. 

(3) Let VEV's of the respective neutral components of = 1, 2) be (l/y/2)pi(z)e l6i( - z \ 
The two moduli pj(z)'s are assumed to take the kink shape of a common width 1/a: 

Pi(z) = v cos/3(l + tanh(az))/2, Pi{z) = t>sin/3(l + tanh(a,2))/2, (2.1) 

where v cos (3 and v sin (3 are VEV's of $1 and $2 respectively in the broken phase 
limit. 

Here the same convention as I is used for the parameters in the effective potential. The 
condition (1) may not be specific since 5 is induced from the soft-SUSY-breaking param- 
eters in the MSSM0. 

Following I, we postulate the effective potential K//, which is considered to include 
the radiative and finite-temperature corrections, as follows: 

11 A A 

V e ff(pi,p 2 ,9) = -mlpl + -mlpl + mlp 1 p 2 cos(5 + 6) + —pl + YPi 

^3 — A 4 2 2 A 5 2 2 Q/i l/\ 2 , \ 2\ a 

+ — ^ P1P2 ~ yPiP 2 cos2# - -{Xep! + \ 7 p 2 )pip2Cos6 

Ap\ + Bp\p 2 cos 9 + Cp x p\ cos 9 + Dpfj , (2.2) 

where 9 = 9\—9 2 - Here the p 3 terms just above are expected to arise at finite temperatures 
so that the kink shape moduli fl2.1|) of the bubble wall are realized for 9(z) = and 5 = 
at Tc- Then several relations among the parameters in (|2.2|) are required as given in I. 

In terms of dimensionless coordinate y = (1 — tanh(az))/2, the equation of motion for 
9{y) derived from V e ff(pi = v cos/3(l — y), p 2 = v sin/5(l — y), 9{y)) is: 

y ( l -y) A 2 + y^ 1 -v)0-- A y)—, — 

dy ay 
= 6sin(5 + 9{y)) + [c(l - yf - e(l - y)\ sin 9{y) + |(1 - yf sm(29(y)), (2.3) 
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where 

b = 



ml 



c = 



4a 2 sin (3 cos (3 ' 

(Aj cot 2 /3 + A 2 tan 2 p + 2(A 3 - A 4 - A 5 )) - — -J= =-r 

2 sin p cos 2 P 



v 2 



32a 2 

,2 



V* 
8^ 



d = 



(A 6 cot/? + A 7 tan/3), (2.4) 
4a 2 ' 

e = . 9 . 2 a ^\Acos 3 p + Dsm 3 p-—) 

4er sm p cos^ p \ v J 

/_B_ g \ 

4a 2 y sin P cos P J 

In addition, the requirement that (pi, p 2 ) — (0,0) and (pi, p 2 ) — (^ cos p 1 , v sin /3) to be 
local minima of V e // with 6 1 = leads to inequalities among the parameters for 5 = 0: 

6>-l, 6 - 2e + 3c> -1 + (A 3 - A 4 - A 5 )^ 2 /4a 2 . (2.5) 

Note that the explicit CP violation 5^0 breaks the symmetry 9{y) < — > —®{y) of 



(P-3|), which is allowed in the case of 5 = 0. 



3 Boundary Conditions Satisfied by 6 

Broken phase limit 

Suppose that, at y ~ 0, is given as 

9(y) = 6 + a y" + (h.o.t.(y)) > 0, a ^ 0), (3.1) 

where (h.o.t.(y)) means (higher order terms of y). Inserting this into fl2.3|) , we have 

y> 2 a + (h.o.t.(y))] = [tyo + W^ + WV] 

+ j,"[W 3 ao + (h.o.t.(y))] (3.2) 
+ ^[(-a 2 /2!)^ 4 + (h.o.t.( 2/ ))], 

where the y v terms on the right hand side come from sin(0(y) — 6q) an d the y 2v terms 
from cos(#(y) — 6> ) — 1- 
That v > requires 

W = b sin(<5 + ) + (c - e + d cos O ) sin 6> = 0, (3.3) 
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from which 

bsind . 

^e = — — — -. (3.4) 

o cos o + c — e + a cos 9 

Without loss of generality, let us restrict 9q as —n/2 < 9 < 3ti/2. Because of b oc m 2 ^ 



necessary to introduce 5, 9q ^ 0, n. Making use of ( |3.3|) 



Wi = (-2c + e-2dcos9 )sm9 = 2bsm(5 + 9 ) -esinOo, (3.5) 

W 2 = (c + dcos^o) sin# = — osin(<5 + 9 ) + esin# , (3.6) 
W 3 = b cos(5 + 9 ) + (c - e) cos 9 + d cos(26> ) 

= — <i sin 2 6*o — &(sin5/sin#o), (3.7) 
W4 = osin(<5 + O ) + ( c ~ e ) sin 6*o + 4d sin #0 cos #0 

= Sdsin^ocos^Q. (3.8) 

If we take v = 1 for illustration, ao is determined as 

q _ 2osin(o" + gg) - e sin go g ^ 
1 + d sin 2 6*o + 6(sin 5/ sin O ) 

from z/ 2 ao = Wi + W 3 a . When v is not an integer, ao is not determined from the lower 
order relations. 
Symmetric phase limit 

Suppose that, at ( = 1 — y ~ 0, 9(Q is given as 

0(C) = 0i + &oC + (h-o.t.(O) (ji > 0, 6 ^ 0). (3.10) 
Inserting this into (|2.3j), we have 



C^(^ + 2)6 + (h.o.t.(C))] = [U + C/iC + ^C 2 ] 

+ CWo + (h.o.t.(C))] (3.11) 
+ C 2At [(-^/2!)t/ 4 + (h.o.t.(C))]. 

That fi > requires £/ = 6sin(o" + X ) = 0, so that 

1 = eic-6 (£ = 0,±1,±2,- • ••). (3.12) 

Making use of this, 

£7i = -esin^ = (-lfesino", (3.13) 

C/ 2 = (c + c/cos^i) sin 6»i = -((-l) e c + dcosd) sin 5, (3.14) 

U 3 = bcos(5 + 9 1 ) = (3.15) 

U A = U = 0. (3.16) 



If we take fi = 1 for illustration, 60 is given from + 2)bo = U\ + C/360 as follows: 

(i) For e ^ and b ^ (-1)^3, 

_ (-l) l esm8 
bo ~ 3 - (-1)^6 ■ (3 ' 17) 

(ii) For e = 0, 6 = (—1)^3 while 60 is not specified apriori. 

Our task is to solve the nonlinear and inhomogeneous differential equation (|2.3j ) with 
the boundary conditions of #0 and 61 (the two-point boundary value problem). We take 
the simplest boundary conditions that 9 ~ 0(5) and 6\ = — 5 throughout, which will 
give the lowest energy configurations. 

4 Examples of Solutions 6 ~ 0(5) 

Assuming that 8(y) ~ 0(5) in the interval y G [0, 1], let us linearize (|2.3|) as (§] 



y\l- y y—^ + y(l-y)(l-±y) 



dy 2 " dy 
= b+(c + d)(l-y) 2 - e(l - y)\e{y) + b5. (4.1) 

Taking v = 1 for illustration, suppose, for a moment, to search for the solution as the 
initial value problem: 9(y) starting from 6(0) = 6$ with 6'(0) = has to hit against 
9(1) = 6\. A naive guess would be that the solution 6(y) should be of 0(5), since 6q 
by fl3.4p, a by ( |3.9p and 6i = —5 are all of 0(5) in general. Actually, numbers of 
such solutions to the linearized equation are obtained even analytically. We show a few 
examples. 

Example 1. We assume 9(y) = po +piy. Given (b, d), a set of algebraic linear equations 
for (po,Pi) obtained from ([4.1|) are satisfied for (c = 4 — d, e = (7 ± y/1 + 166 2 )/2), giving 
Po = ~(b/(b + c + d — e)) x 5 and pi = — (4/(e — 3))po- F° r this 9(y) to satisfy the 
boundary condition 9 1 = —5, the unique choice is (b,e) = (3,0). Namely, for (b,c,d,e) = 
(3, 4 — d,d,0), we have a linear solution 

9(y) = -((3 + 4y)/7)x5. (4.2) 

Example 2. In the similar way, for (b,c,d,e) = (3,10 — d, d, 0), we have a quadratic 
solution 

%) = -((3 + 5y + 5y 2 )/13)x5. (4.3) 

In the case of spontaneous CP violation, we obtained an almost linear solution 9 s=0 (y) 
shown in Fig.l of I, for which —9 s=0 (y) is also a solution. However, once 5 ^ is given, 
—9(y) of ( |4.2[ ) or ( |4.3| ) is no more a solution. Although we have no cancellation in the 



net baryon number in this case, such solutions of 9(y) ~ 0(5) would contribute at most 
marginally to the baryon asymmetry f|]]. 
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Figure 1: Numerical solutions of 9 5=0 (y) for (&(°), c (°\ S°\ e<®) = (3, 12.2, -2, 12.2) (solid 
curve), (3,8.98,1,11.98) (dashed curve) and (3,10,-0.2,11,8) (dotted curve). The first 
is the one given in Fig.3 of I. For these, 9o = 9± = 0. 

5 Solutions 6 ~ 0(1) around Bubble Wall 

Suppose that for 5 = we have a set of solutions ±9 s=0 (y) which are not 0(5) but 0(1). 
For 1 5 1 small enough, we could expect a solution 9 + (y) close to 9 5=0 (y) and another one 
0~(y) close to —9 5=0 (y). The both 9 ± (y) satisfy the same boundary conditions 9 ~ 0(5) 
and 9i = —5 but 9~(y) ^ —9 + (y). Before giving numerical solutions 9 ± (y), we show 
how such solutions of 0(1) are obtained from ±9 5=0 (y) for sufficiently small \5\. For 
definiteness, we take v — /j, — 1. 
Solutions ±9 5=0 (y) 

In Fig.l we show three 9 5=0 (y) satisfying the boundary conditions 9 = 9 X = together 
with the corresponding parameters (b^°\c^°\S°\e^). One of them is the one given in 
Fig.3 of I. Because W± — U\ — for 5 — 0, the parameters are required to satisfy 

6 (0) + c (0) _ e (0) + d (0) = u 2 = ^ (51) 
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from [y 2 - W 3 \s=o)a 5 ~° = 0, and 

&(°) = (ifa + 2) = 3 (5.2) 



from + 2) — U 3 \s=o)bo = 0. Of course, once we get a solution 9 s °(y), its a s 
and &o = ° have been uniquely determined to match the boundary conditions. As is clear 
from the figure, 0'(O) = a 5 =0 ~ 0(1) and 0'(1) = -b 5 =0 ~ 0(1). These are what enable 
9(y ~ 0.5) to become as large as 0(1). 
Solutions 9 ± (y) 

This suggests that we could have solutions 9 ± {y) ~ 0(1) for 5 ^ if a ~ 6 ~ 0(1). As 
shown in Section 3, they are determined from the lower order relations when v — \i — 1. 
Now we seek a set of parameters (b, c, d, e) , which incorporates such solutions for v = // = 
1.0 

Put 6 = + Abx5,c = c (0) + Ac x 5, d = d (0) and e = e (0) + Ae x 5. For sufficiently 
small 1 5 1, we have from (|3.4|), with the use of (|5.1|) , 

O ~ -i : x 5 ~ -6 x 5, i.e., 5/9 ~ -1/6. (5.3) 

This implies that the denominator of a® in ( |3.9| ) is not of 0(1) but 0(5). Since the 
numerator of clq is 0(5), ao is now a quantity of 0(1): a ~ 6^°) (2c < - - ) — e^ -* +2<i^)/(— A6 — 
Ac + Ae). Because of ( |5.2|) , 6 in ( |3.17| ) is also a quantity of 0(1) for e^ ^ 0: 6 ~ 



—e^/Ab. By suitably adjusting (A6, Ac, Ae) in such a way that ao ~ ag ~ 0(1) and 

&o ~ &o = ° ~ 0(1) match the boundary conditions, we obtain a desired solution 9 + (y) close 

to 9 s=0 {y). For the same ( A6, Ac, Ae) or (6, c, d, e) , we can find another desired solution 

9~{y) close to —9 5=0 (y). In Fig. 2 we show an example of 9 ± {y). Note that —9~(y) 

and 9 + (y) do not coincide with but considerably differ from each other. Note also that, 

while 9 ~ — b x 5 and 6\ = —5, 9 ± (y) deviate nonperturbatively from the corresponding 

±9 s=0 (y) in the intermediate region. Presumably this may be due to the nonlinearity and 

the singular effects for 9"{y) near y = 0, 1 of the differential equation (|2.3|) . 

We can also find several solutions for other sets of the parameters (b,c,d,e), as long 

as they do not change the global structure of the effective potential. Of course, they 

would not have v = fi = 1 in general. The numerical method is based on the relaxation 

algorithm. For sufficiently small S, say 5 G (0, 0.001], the both types of solutions 9 ± {y) can 

be found by starting from initial configurations with opposite signs, respectively. As we 

increase S, only 9~{y) can be obtained starting from any initial configuration. This would 

2 In practice, we are given with the effective potential, so that (b, c, d, e) are fixed. Then v and [i are 
determined by solving the equation for 9(y). Since our purpose here is to show the possibility to have 
such solutions, we trace an unusual course to find them. 
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Figure 2: Numerical solutions of ± {y) in which 5 > is incorporated into the symmetric 
pairs of the solid-curve solution 9 6=0 {y) for (&(°), c<~°\ S°\ e<®) = (3.12.2, -2, 12.2) in Fig. 1 
by the prescription explained in the text. These pairs are given respectively by the upper 
and lower solid curves. For 5 = 0.0025 and (b,c,d,e) = (2.98005,12.178375,-2,12.2), 
0~(y) is given by the dashed curve and 6 + {y) by the dotted one. They have the common 
boundary values: 0jj= = -3.109066 x 5 x 5 and 9 1 = -5. 
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mean that for larger 5, 9 + {y) becomes more unstable or even does not exist Obviously 
the critical value of 5, above which 6 + (y) is not found, depends on the parameters in the 
potential (b, c, d, e). Although its value is about O(10~ 3 ) ~ O(10~ 2 ) in all the cases we 
studied, we could not determine its definite value, since it might depend on details in the 
algorithm, such as the convergence parameter or initial configurations. 



6 Energy Density and Enhancement Factor 

The energy density of the wall per unit area is given by 



S 



oo I 1 

* i=l,2 



' dpi 



dz j ^ l l dz , 



dy { ay(l -y) 



1=1,2 



dp 



+ K//(Pl,P2, 
2" 



(6.1) 



dd,; 



-r 1 +Pt -n 



dy 



dy 



+ 



203/(1 - y) 



V eff (p 1 ,p 2 ,0) 



= av 2 /3 + £[9], 



(6.2) 



where the first term of the above line is the energy density of the trivial solution 9 5=0 (y) = 
0. The second term is contributed from a path 9{y) that connects the broken and sym- 
metric vacua. For the trivial solution this term vanishes. 

£ has two degenerate minima corresponding to ±9 5=0 (y) ^ 0. For the solid-curve 
solutions in Figs.l and 2, we have 



s[±e 



<5=0i 



-2.056 x 10-Wsin 2 /3cos^ 



(6.3) 



That (5/0 breaks the degeneracy. For 6 ± (y) shown in |2|, in which 5 = 0.0025 is 



incorporated into these ±6 5 -°(y), we actually find that > S[6 + ] > S[±6 s -°] > S[6 
the energy difference being 

AS = £[9~\ - £[9 + ] = 



1.917 x lO'Vsin^cos*/?. 



(6.4) 



This negative A£ means that the formation of the bubble with 9 (y) is favored over that 
with 9 + (y). The relative enhancement factor is given by 



exp 



AttR 2 c A£^ 



T, 



c 



(6.5) 



where the radius of the critical bubble Rc is approximately given by \j3Fc/(4:7rav 2 ) with 

Fq being the free energy of the critical bubble. Various authors estimate Fc ~ (145 ~ 

3 This disparity in 8 ± {y) would be triggered by #o,#i < for 5 > in our choice and amplified 
nonperturbatively as remarked above. 
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160)T. If we take Fc = 145T and tan/3 = 1, the enhancement factor is 

exp jg- = 8.05. (6.6) 



c 



Such a large relative enhancement factor of 0(10) between the bubble with 9~(y) and that 
with 6 + (y) would surely avoid the cancellation in the chiral charge flux and guarantee the 
survival of the net baryon number of the universe. 



7 Concluding Remarks 

Of course, the simplified condition (3) at the beginning of Section 2 that the wall moduli 
are fixed to be the kink shape is invalid for the solutions 9(y) ~ 0(1), and we have to solve 
a set of coupled equations for 9{y) and Pi(y) as done in [[J. As shown there, 9{y) remains 
to be 0(1) though its form is modified to certain extent while Pi(y) 7 s largely deviate from 
the kink shapefj The kink-shape approximation is valid only for the solutions of 0(5) 
as given in Section 4. Though they have no counter partners as remarked there, the 
net baryon number would remain at most marginally after the completion of the phase 
transition. 

On the other hand, as we gave an estimate in I for the solid-curve solution 9 s=0 (y) in 
Fig.l, such 9 ± {y) that is able to become 0(1) around the bubble wall could supply an 
efficient chiral charge flux through the wall surface at the phase transition. Note that in 
the broken phase limit the CP violation is given by sin(5 + 9q) ~ 5 at Tc, so that, if \5\ 
is small enough, there should be no contradictions with the present experimental bounds. 
Of course the CP violation completely vanishes in the symmetric phase limit because of 
Pi = there. That solutions with such features as presented here are allowed in a realistic 
model would be highly significant in any scenario of the electroweak baryogenesis, since 
there may be numbers of possible mechanisms to diminish the net chiral charge or the 
net baryon number before the completion of the phase transition. 

Our estimate in Section 6 suggests an interesting possibility that the effect of a small 

\S\is nonperturbatively amplified to yield a large relative enhancement factor to favor the 

formation of only one of the two kinds of the bubbles, which are the symmetric partners 

in the absence of the explicit CP violation. Further our numerical analysis suggests that 

the bubble with higher energy would be metastable or could not exist for larger \S\. If 

the possibility is realized, we could be free from the disgusting cancellation in the baryon 

number of the universe. 

4 For a numerical method how to obtain the chiral transmission and reflection coefficients in such cases, 
see §. 
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